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DEVELOIMENT OF A LAMINAR BOUNDARY LAYER 
BEHIND A SUCTION POINT* 

By W. Wuest 


1. INTRODUCTION 


Boundaiy- layer suction originally was applied to reduce the 
boundary- layer thickness and therewith the inclination to flow sepa- 
ration; however^ since the properties of bodies with small drag have 
been improved more and more, attention was drawn to an increased extent 
to the reduction of surface friction. One now strived toward keeping 
the boundary layer laminar as long as possible, thus to defer the tran- 
sition point to turbulence as far as possible. Boundary- layer suction 
was recognized to have a favorable effect in this sense, and therewith 
the velocity distribution in a laminar boundary layer behind a suction 
point acq,uired hei^tened interest. The stability of a laminar velocity 
profile is very severely affected 'by the shape of this profile. 

In a considerable number of theoretical reports (reference l) the 
case of continuous suction was treated for reasons of mathematical 
simplicity; permeability of the wall surface was assumed. In further 
reports, the stability of laminar boundary- layer profiles in case of 
continuous suction was treated and a considerable rise in the stability 
limit was determined; however, a technical realization of such peme- 
able walls with sufficiently smooth surface and adeq.uate material 
strength characteristics is difficult. For structural reasons, it is 
simpler to arrange single-suction slots. In addition to the suction 
effect proper, there appears here the sink effect first discussed in 
detail by L. Prandtl and 0. Schrehk (reference 2 ) and recently treated 
by Pfenniger (reference 3) in an instructive experimental Investigation. 

Below, the pressure variation along the wall as well as, in partic- 
ular, the sink effect are disregarded. Figure 1 shows the practical 
realization of such a case. We assume that on a flat plate A, a laminar 
boundary layer ("Blasius boundary layer") develops at constant pressure. 
We assume a second plate B arranged beginning from a certain point Xq 

at the distance y^ parallel to the first plate so that a suction slot 


"Entwlcklung einer laminaren Grenzschicht hinter elner 
Absaugeste3JLe.” Ingenieur Archiv, Vol. 17^ 19^9^ PP* 199-206. 
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is formed between the tvo plates. The magnitude of the power req.uire- 
ment for suction is assumed to be precisel;^ such that^nerely the part 
of the boundary layer situated between the ^two plates is removed. Thus, 
there begins above the plate B a new laminar boundary _^layer which is. 
distinguished from the Blaalus boundary layer by another Initial condi- 
tion. The new boundary layer foims at its ^tart the .outer part of a 
Blaslus boundary layer. 


2. BOUMDARY-LAYEE EQUATION AND ASYMPTOTIC BEHAVIOR 


By introduction of the stream function_and the total pressure, the 
boundary- layer eq,uatlon may be transformed by the well-known method 
(reference 4) into 


5x 


VU 




( 1 ) 


where g = p + £ u^ and ^ = u. We limit ourselves to the case that 

oy r 

the flow takes place outside of the boundary layer at the velocity 
U]^ = const. , thus to the flat plate and put furthermore 


g = - | + Const. 


or, respectively 


u = 


( 2 ) 


(3) 


This statement has been chosen so that for large ij;-values, q assumes 
the value 1. Equation (l) is thereby transformed into - 






( 4 ) 


From the definition of the stream function and from equation 


further obtains with 



(3), 


one 
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- (5) 
. ( 6 ) 


y = 



(7) 


In order to investigate tlie asyn^jtotic behavior of the differential 
eq.uation (3)^ we put for large values of -^r 


q. = 1 “ q.-vr with 0:^^ « 1 


(8) 


In first approximation, one then obtains 

8q, 


Sx 


^ = vu.2_ 




(9) 


This differential eq.uation, however, is mathematically identical with 
the differential equation of a nonsteady flow independent of x which 
has been treated before (reference 5)^ the time t is now replaced by 
the stipulated space coordinate x. It also corresponds to the well- 
known heat-conduction equation. The general solution is therefore given 
by 


5 ,/ t + r \ 


■ ^o)j 

- *o)y 


r q^(0,x»)^ $f-== 

^ xq 8x ^^kvu^ 


t > 

(X - X’); 


\dx‘ 


dilr’ + 


( 10 ) 
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Therein 


0 = r e-y^ dy 

is the kncfwn error Integral. W. Tollmien (reference Q. has investigated 
this solution for two special cases where the first integral disappears. 
For the boundary layer with suction, however, this will no longer be the 
case. 


3.. BLASIUS BOUNDARY LAYER 


Althou^ we presupposed that the velocity U]_ at the edge of the 

boundary layer is constant, the problem of the suction boundary layer 
to be treated here nevertheless differs from the flow on a simple flat 
plate ("Blaslus boundaiy layer”) by the fact that othe:^lnitial condi- 
tions exist; rather, the Blasius boundaiy lajrer is contained as special 
solution among the suction boundary layers since there "Xq = 0, thus 

suction point and beginning of the plate A (fig. l) coincide. Since we » 

shall make use of this special solution for the later calculation, we 
shall first consider the Blasius boundary layer. It is _diBtinguished ^ 

by the fact that q may be regarded as dependent merely on a quan- 
tity q = One then obtains from equation {h) the following 

differential equation of the Blasius boundary layer ~ 






( 11 ) 


The solution may be written in the following form 


q.B = PTI 




( 12 ) 
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The constants therein have the following values 



ao = ! 

^ (990 X 15) 

ajj^ = 1.60333 X 10 “^ 


a^ = 0.57627 X 10 "^ 

a^ = 3.8907 x'l0“9 

A 

ay = -1.3986 x ' 10-9 

ag = -3.9135 X 10"^^ 


a^ = 3.7282 X 10 "^^ 

= 2.2383 X 10’ ^3 
an = -0.3101^ X 10-1^ 
a^2 = -1.081 X 10-15 


Due to the boundary condition at tbs wall, one integration constant is 
zero. The second integration coi;^ 8 ta^t is determined frcm the asymptotic 
behavior for large values of t]. Because of q^Cijr^O) = 0 the first 

Integral in eq.uation (lO) is elimi'f^ted. The second integral, however, 
yields by partial integration, with consideration of the asymptotic 
behavior of the error Integral, just as in the case treated before by 
W. Wuest the solution 





(13) 


The constant^ P in eq.uation ( 12 ) and 7 in equation (13) are deter- 
mined by the fact that for large -q values q and Sq/dTj according 
to equation (12) and equation (I3) agree ^ith eaph other. Th^ recalcu- 
lation of the two const^ants yleld(^ the following values 

P =^ 0 . 6642 ' . .7 =^0.828 

\ 

For comparison, L. Prandtl (reference 7 ) gives^'thj f6ilowing values 
calculated by Blasiqs an^.Tg^Ufer which read, converted to the above 
designations . ' jk ' 


P = 2 X 0.332 = 0.664 

7 = 2/jT 0.231 = 0.819 
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F. Elegels and J. A. Zaat give in a new report (reference 8) for 7 
the following value — 


7 = 0 . 3 ^ 2 \^ = 0.857 


The function q with first and second derivative has been tabulated 
and plotted in numerical table 1 and figure 2^. 

Ij-. ASmPTOTIC BEHAVIOR OF THE SUCTIOH BODNMRY lAYER 


For calculation of the asymptotic behavior of the suction boundary 
layer, we divide the function defined by equation (8) into two 

parts 

“ S^l Sf2 


The first part is to be selected so that it satisfies the initial condi- 
tion at the suction point x = Xq-j this is done.' by extending the asymp- 
totic solution of the Blasius boundary layer to. x > xq as well. From 
equation (I3) one then obtains 


^1 = ^ 


1 - 3>i 


+ VqN 


and 


^he numerical table has been calculated with the values 

7 = 0;819. 


3 = 0.664 
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Numerical Table 1. Blasius Boundary liayer 


n 


q.^(Tl) 

d“(n) 

\/ci(n) 

0 

0 

0. 664o 

0 

0 

.1 

.06606 

.6555 

.12751 

.2570 

.2 

.13106 

.6427 

.17750 

.3620 

.3 

.1939 

.6206 

.2114 

.4404 


.2546 

.5981 

.2369 

.5046 

.5 

.3135 

.5734 

.2561 

.5599 

. 6 

.3695 

.5471 

.2700 

.6079 

.7 

.4228 

.5195 

.2796 

.6502 

.8 

.4681 

.4912 

.2856 

.6842 

.9 

.5211 

.4625 

.2883 

.7212 

1.0 

.5659 

.4337 

.2883 

.7523 

1.1 

.6081 

.4062 

.2857 

.7798 

1.2 

.6469 

.3766 

.2809 

.8043 

1.3 

.6831 

.3487 

.2743 

.8265 

l.k 

.7167 

.3217 

.2660 

.8466 

1.5 

,7474 

.2955 

.2563 

.8645 

1.6 

.7758 

.2705 

.2457 

.8807 

l.T 

.8017 

.2466 

.2340 

.8954 

1.8 

.8252 

.2238 

.2248 

.9084 

1.9 

.8465 

. 2022 

.2088 

.9200 

2.0 

.8657 

.1820 

.1955 

.9304 

2.5 

.9352 

.1013 

.1303 

.9671 

3.0 , 

.9715 

.0509 

.0774 

.9857 

3.5 

.9881 

.0217 

.0381 

.9940 

4.0 

.9962 

.0085 

.0169 

.9981 

4.5 

.9988 

.0029 

.0066 

.9994 

5.0 

.9997 

.0009 

.0022 

.9998 


Therein, ijr = 0 forms the nev wall streamline and tItq the suction 
q.uantity. The secoQd part q ^2 then must be chosen so that the bound- 
ary condition = q^(0,x) is satisfied. If the asymptotic rela- 
tion q. ~ 1 - q.^ would rigorously apply in the entire domain of the 
boundary layer, there would have to be at the wall qy(0,x) = 1, because 

of q. = Oj however, the asymptotic soluticai deviates from the rigorous 
solution if it is continued up to the wall. Therefore q^(0,x) is an 

untoiown function regarding which we merely make the assumption that it 
does not become infinite. As initial condition for the. part q:^^ one 
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further has ^ ° since q^i(t,XQ) already satisfies the 

initial condition 


q-v = 7 


n , ’■I'n 

1 ^ 0^ 


Iv^i^o/ 


which insures the connection w^th the Blasius soluticai. The contrl- 
hution to thi^ solution al-afb must obey the differential equa- 

tion (9). In the solution (equation (lO)} the first integral is ellml- 
nated, because of whereas in the second Integral one 

has to put • ■ 


= q^(0,x) - q^3_(0,x) = q^(0,x) - 7 
so that the asymptotic solution reads 


1 - $/ 




/hvu^x^ 


= r 


1-4 




- r 


1 - 




yhvuix' 


5*’ t^4vui(x - 


By partial Integration one obtains with consideration of the asymptotic 
behavior of the error integral (by W. Wuest, else-vdiere) 


1 - $/ 


'If + \ 
l^ 4 vu^ 


\~l 




- 7 


1 - ® 




1 - (t/ 


|/ 4 VUi(x - Xq) 


Because of the connection with the Blasius solution, however, 
q^(OjXo) =7, if the asymptotic solution is continued up to the wall, 

so that one finally obtains as the asymptotic, solution for the suction 
boundaiy layer 


Sf ~ 7 



+ 7i>f \ 

1 - o/ ^ 

\v^4vuixy 

V4vuiXo/ 

1 

M 

1 

' ! 


(14) 
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Instead of the error Integrals $ one may for large values of "'If 
again go back to the Blasius solution if one takes the asymptotic 
behavior of the latter according to equation ( 8 ) and eq.uatlon (I3) into 
consideration 

r-g j , 

yvui(x - xq ) 

In this formula represents the Blasius solution. The last form of 

the solution proves to be particularly expedient for the further con- 
siderations. 




+ ^0^ 


- -5 


W^VuiXoj 



5. APPROXIMATE SOLUTION FOE THE SUCTION BOUNDARY LAYER 


It suggests itself to generalize the asymptotic solution which is 
valid for. large values of in the following manner 


q = 



F(ilf,x) 


1 - q 






/VU1(X - Xq)^ 


(16) 


Due to q = 0 for -ijf = 0 and because of equation (15) the func- 
tion F(\!f,x) must satisfy the following conditions 


F( 0 ,x) 



F(oo,x) 




fhVuiXo 


(IT) 


It was hoped at first that one could choose for F, as in the nonsteady 
analogue by W. Wuest, elsewhere correspondingly an exponential func- 
tion as the simplest formulation; besides equation (IT) th 5 ~dls appear- 
ance of the second derivative of q at the wall would be added as a 
further condition; however, it was shown that such a formulation does 
not meet with success and even, in a certain domain, does not yield any 
solution at all. 
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For the further calculation we Introduce the following simplified 
notation 


l/VU3^(x - Xq) 

BO that the solution (equation 

<1 = q3(n) - rjjl - q^(il')jj (19) 

According to a suggestion by A. Betz, we equate as first approxi- 
mation q^ the function F to the value dependent only on x 


_Jc 


ti = no + 






n’ (18) 


(16)) reads 


Fq(x) 


F(0,x) 



= ‘iB(no) 


at the wall. * Thus the first approximation reads 


(20) 


q^ = qg(ri) - Fq^I - Q-gCn’)^ (21) 

This formulation does not fulfill the boundary- layer equation (it) exactly. 
In particular, the second derivative of q^^ "at the wall does not disap- 
pear; however, the dependency on the second derivative of the stability 
of the velocity profile Is of a very sensitive nature, eo that one has to 
look for a more accurate solution. By substitution of the approximate 
solution '(equation (21)) Into the boundary- layer equation (4), one 
obtains 


V‘U)- 


- ^ 

X 




Hence there results as the error of this first approxi- 

mation. By subtraction of the exact solution .In which F stands 
for Fq and q for q^, while disappears, one then obtains 
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(22) 


■Mhere 




;i = . 2 " . 2(- - =o)a^ 




,, a. uBkno» function. Tha auaatlty V' ° 

. V . „. By integration of elation (22) one obtains 


Therein S3_ is i*® 
quadrature 


(P - Po)[^ ■ ®2 

detemlned graphically or numerically hy repeated 


pTl’ p '^2 ^^61 A 

- 'I L ^ 




(21^) 


One may determine the asymptotic behavior of 

above definition of Sg*' for ^ end \/J 

values 


gg hy substituting in the 
the asymptotic 
Thereby one obtains 


gg« « ~ (x - 


Bence foUoVB vlth use of equations (9), Ug). (21), 
gration vltb respect to n' = 


and repeated inte- 
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Aa before, $ denotes the error integral. Generally we visualize ^2 
as represented in the following manner 


^2 = 21 





( 26 ) 


By way of approx^Lmation we limit ourselves to the first two terms, 
with a-j_ = - Fq) and a^ determined hy the fact that q. must 

disappear at the wall. We determine accordingly the faction F 
approximately to be 


F = Fq + 


1 - 


^bU') 


. 6 ^(n',x) + r(F„ - Fo)|i - + a^ll - ®‘(n‘3l| 


(27) 


a^ = -6^(0, x) - r(F<„ - Fq) 


( 28 ) 


Calculation example .- = 0,125 was selected as numerical 

exangjlej F was calculated for the values x/xq = 1.234, I. 562 , 4.34, 
and 9«?8 and plotted in figure 3 . For x/xq = I .562 the error was 

determined by substitution of the approximated solution .into the 
boundary- layer equation, and compared with the- first approximation 
according to equation (21). Con^iared to the first approximation, a con- 
siderable improvement results particularly in the region^near the wall 
(fig. 4). In figure 5 "the results are converted to the velocity pro- 
file, in figure 6 the second derivative is represented. As a supplement, 
the connection between the degree of suction and the suction quantity of 

the magnitude tJq* = will be supplemented. By the degree of 

suction 0 we here understand 


0=1--^ (29) 

B]_* being the displacement thickness immediately ahead of the suction 
point and & 2 * immediately behind it. (Therewith, 0 is given by 
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(30) 


The resulting values are tabulated in table 2 and plotted in figure 7 * 


Table 2 . Degree of Suction 


a 

0 

0.1 

0.2 

f 

0.4 

0.6 

0.8 

1.0 

2.0 

© 

0 

0.392 

0.520 

0.671 

0.762 

0.82k 

0.870 

0.97^ 


The suction q.uantlty i|fQ Is^ furthermore, given by the following 
relation 

(31) 

6 . SUMMARY 

The developmen^of st laminar boundary layer behind a suction point 
is investigated if b^the suction merely the part of the boundary layer 
near the wall is "cut off", withcjut the slot exerting a sink effect. 

As basis of the calculation, we used the boundary- layer equation in the 
form indicated by Piundtl-Mlses which is . closely related to the heat 
conduction equation oh, respectively, to the differential equation of 
the nonsteady flowCHdilch Is^ independent of the coordinate x along the 
wall. With consideration of '^h^asympt'otic behavior of the solution, 
an approximate solution is developed which is similar in structure to 
the solution of the nonsteady analogue which has been treated in an 
earlier report by W. Wuest, . elsewhere. 


Translated by Mary -L, ^ahler 
Rational Advisory Conmlttee 
for Aeronautics ■ 


.1 
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Suction 


Figure 1 .- Boundary -layer suction at the flat plate without sink effect. 



Figure 2 ,- The function q(T]) of the Blasius boundary layer with first and 

second derivative. 
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Figures.- Axixiliary function F(n*,x) for calculation of the suction 
boundary layer for ^Q/yvn^XQ = 0,12fc. 



Figure 4.- Error of the first and second approximation for = x/xQ = 1,562. 
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Figure 5.- Velocity profiles of the suction boundary layer for 
" 0*125 ^d various distances from the suction 

point. 



Figure 6,- Second derivative of the velocity profiles of the suction 
boiindary layer for = 0.125. 
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Figure 7.- 



Degree of suction (ratio of the cross-hatched and the total 
shaded area in fig. 1). 
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